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Abstract Supergranulation is one of the most visible length scales of solar 
convection and has been studied extensively by local helioscismology. We use 
synthetic data computed with the Seismic Propagation through Active Regions 
and Convection (SPARC) code to test regularized-least squares (RLS) inversions 
of helioseismic holography measurements for a supergranulation-like flow. The 
code simulates the acoustic wavefield by solving the linearized three-dimensional 
Euler equations in Cartesian geometry. We model a single supergranulation cell 
with a simple, axisymmetric, mass-conserving flow. 

The use of simulated data provides an opportunity for direct evaluation of the 
accuracy of measurement and inversion techniques. The RLS technique applied 
to helioseismic-holography measurements is generally successful in reproducing 
the structure of the horizontal flow field of the model supergranule cell. The 
errors are significant in horizontal-flow inversions near the top and bottom of 
the computational domain as well as in vertical-flow inversions throughout the 
domain. We show that the errors in the vertical velocity are due largely to cross 
talk from the horizontal velocity. 

Keywords: Interior, Convection Zone; Supergranulation; Helioseismology, In- 
verse Modeling; Helioseismology, Direct Modeling; Waves, Propagation; Velocity 
Fields, Interior 



1. Introduction 

Wave travel times are influenced by the properties of the medium in which 



they propagate; helioseismic holography (Lindsey and Braun 1990 2000) can 



be used to measure travel-time shifts, and is therefore useful for probing the 
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interior structure of the Sun ( Gizon and Birch 2005 Gizon, Birch, and Spruit 



2010) 



Supergranulation is a visible pattern on the solar surface often interpreted 
as a scale of solar convection. Despite extensive research by local helioseismol- 
ogy, there remains substantial uncertainty regarding the characteristic scales 
and fundamental processes responsible for the creation and sustainment of su- 
pergranulation patterns (e. 



2010). 



Gizon and Birch 2005 Rieutord and Rincon 



Supergranulation cells were first detected by Hart (1954 1956); observations 
and analysis of data over a period of five decades have revealed average cell 
sizes in the range 15 - 30 Mm (e.g. Hart 1956 Leighton, Noyes, and Simon 



196 2] |Duvall[ [1980) |Ha~thaway[ [19921 |Hagenaar, Schrijver, and Title! |1997[ |De 



Rosa and Toomre] |2004| |Hirzberger et aL||2008 ), with discrepancies in estimated 
scale due largely to differences in analysis techniques. Similarly, direct inferences 
of the velocity field have produced a range of horizontal-flow estimates (e.g. 



Hart 




1954 


al. ( 


2010 


), 



1954 Simon and Leighton 1964 Hathaway et at, 2002). Rieutord et 



al. (2010), motivated by observed scale-dependency in surface kinetic-energy 



distributions, performed a spectral analysis of velocities estimated using granule 
tracking and Doppler measurements. They measured a velocity of 300 ms' 1 at a 
scale of 36 Mm. Estimating vertical velocity has proven more difficult due to the 
noise contribution and low signal strength, particularly near the boundaries of 



super granulation cells ([Rieutord and Rincon 2010[ ). Results from SOHO/MDI 
data ( |Hathaway et al. |2002[ ) a nd Hinode/ SOT data (|Rieutord et al.\ |2010[ ) 



Duvall and Birch (2010) estimated cell-center 



using direct Doppler measurements averaged 



produced estimates of 30 ms" 
vertical flow to be about 10 m s" 
over 1100 supergranules. 

Several theories have been advanced to explain the development of supergran- 
ulation flow patterns at the observed preferential scale, although confirmation or 
rebuttal of these theories is largely lacking. An early theory proposed by |Simon| 
and Leighton (1964) attributed supergranular scales to convective instability 



from the recombination of ionized helium. More recently, Rieutord et al. ( 2000 



2001 ) proposed a model whereby exploding granules trigger large-scale instability 
of the granular flow, leading to supergranulation. |Ploner, Solanki, and Gadun] 
(2000) used a ID model to similarly show that the interaction and merging 



of individual granular plumes may determine supergranular scales. Rast (2003) 
demonstrated that large spatial and long temporal supergranular scales arise in 
a simplified advective model that simulates the interaction of many small-scale 
and short-lived granular downflow plumes. [Crouch, Charbonneau, and Thibault] 
( 2007 ) were able to simulate correlations between cell size and magnetic activity 



using a model based on a random-walk approximation to the dispersal and 
interaction of small-scale magnetic elements at the solar surface. Qualitatively, 
they found that this process produced supergranule-like spatial patterns. 

Debate continues regarding the structure of supergranule flows below the 



solar surface (e.g. Duvall 1998 Zhao and Kosovichev 2003 Woodard 2007 



Sekii et al. 2007), largely due to the difficulties associated with estimating 
relevant quantities using helioseismology (e.g 



Gizon and Birch 2005 Rieutord 



and Rincon, 2010). These challenges are compounded by the sound speed and 
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density stratification, which effectively decrease resolution and signal-to-noise 
ratio with depth, making the detection of a return flow elusive. |Kosovichev 



and Duvall 



vail et al. 



( 1997 ) first used the time-distance helioseismology technique ( Du 



1993) to make measurements of supergranular flow patterns from 



Doppler measurements. Zhao and Kosovichev (2003) also used time-distance 



helioseismology to infer supergranule flow fields; they found that horizontal 
flows could be derived reliably within a few megameters of the surface, but 
that vertical flows remained uncertain. They determined supergranulation to 
be a relatively deep phenomenon with a full pattern depth of up to 15 Mm 
and a return flow measured below a depth of 5 — 6 Mm. Braun, Birch, and 



Lindsey (2004), using helioseismic holography, attributed the measured return 



flow to signal leakage, and concluded that supergranulation is a relatively shal- 



low phenomenon. Woodard (2007) used a forward model to perform subsurface 
inversions of Doppler data, however detection was limited to 4 - 5 Mm below 



the photosphere Jackiewicz, Gizon, and Birch (2008) used a novel 2+1D opti- 



mally localized averaging (OLA) method to study supergranulation in the top 
few Mm. IBraun et al.\ (I2007I) and IZhao et al.\ (I2007I) used synthetic data from 



numerical simulations of wave propagation through supergranule-like flows to 
test the performance of helioseismic techniques. The former compared model 
travel times with those computed from synthetic observations using helioseismic 
holography. The latter computed inversions using time-distance helioseismology 
and compared the results to the modeled flows. Both studies documented limi- 
tations in the a bility of the techniques to detect the full extent of the flow fields. 



Svanda et al. 



(2011) tested the OLA inversion technique of 



Jackiewicz et al. 



(2012) on synthetic travel times and showed improved ability to infer vertical 



flow fields. 

In this study, we simulate wave propagation through a simple kinematic model 
of a supergranule flow pattern and make helioseismic holography measurements 
from observations of the resulting velocity field. The aim is to test the quality of 



regularized least squares (RLS: e.g. Zhao, Kosovichev, and Duvall 2001 Koso- 



vichev 1996 in the context of local helioseismology) inversions and gain insight 
into the limitations of this commonly used technique in local helioseismology 
by comparing the calculations with a known simulated flow field. We find that 
the RLS technique applied to helioseismic holography measurements is able to 
infer some general features of the horizontal field, but fails to infer the vertical 
field throughout the domain. Errors in vertical-flow inversions have previously 



been attributed to cross-talk effects (e.g. Zhao et al. 2007 Jackiewicz et al. 



2012); herein we show the individual contributions from the divergent flow that 



comprise these effects. 



2. Simulations 

2.1. Numerical Algorithm 



We used the SPARC code (Hanasoge et al. 2008 2006 Hanasoge and Duvall 



2006) to solve the three-dimensional linearized Euler equations in Cartesian ge- 



ometry. The code computes derivatives in the vertical direction using sixth-order 
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compact finite differences (Leie 1992), whiie derivatives in the horizontal direc- 



tions are computed spectrally with periodic boundaries. An absorbing sponge 



is used to damp wave reflection at the top and bottom boundaries (Hanasoge 



et al. 20061. 



Wave excitation is applied 200 km below the photosphere via a source term in 



the momentum equation (Hanasoge, Duvall, and Couvidat 2007). The forcing 



is statistically uniform in the horizontal directions and Gaussian in the vertical 
direction with FWHM 200 km. In the spectral domain, the forcing distribution 
is localized to within the wavenumber range k = — 2 rad Mm -1 and has a 
frequency maximum at 3 mHz. 

Two different 24-hour (solar time) simulation cases were performed. The first 
case used the supergranule model with velocity vector v = v x x + v y y + v z z, 
where x and y define the horizontal coordinate system and z defines vertical 
distance from the photosphere. The second simulation case used no background 
flows, which allowed for application of the "noise subtraction" technique in the 



analysis (Werne, Birch, and Julien 2004) 



2.2. Supergranulation Model 



The static background model is a convectively stabilized (CSMJV: Schunker 
20111 variant of Model S ( Christensen-Dalsgaard et al. 



et al. 



1996) 



which 

specifies the time-invariant components of the density [/?], pressure [p], first 
adiabatic index \Ti], gravity [<?], and sound speed [c s ]. The computational box is 
a 300 3 grid, spanning nonuniformly from Zbot = ^25 Mm to Zt op = 2.5 Mm in the 
vertical direction. The grid spans 100 Mm uniformly in the x- and y-directions. 

We model a single supergranule with a simple axisymmetric mass-conserving 
flow (Figure [IJ prescribed by the curl of a potential function 



pv = V x A, 



(1) 



where, for the sake of simplicity, we assume that the potential is separable as 



A = A^ = P {z)h{z)F{r)cj> 



(2) 



with functions h(z) and F(r) defining the vertical and radial dependance of 
the potential. The variables (f>, r, and z denote the angular, radial, and vertical 
cylindrical coordinates, respectively. The formulation of Equation [T] is chosen so 
that mass conservation V • pv = is automatically enforced. We chose the radial 
distribution F(r) to be of the form 



F(r) = Ji(fcr)e- r/fl , 



(3) 



where J\ is the first order Bessel function of the first kind, k = 27r/30 rad Mm -1 , 
and R = 15 Mm. This formulation forces decaying reversal of flow direction with 
increasing radial distance from the axial center. This functional form is similar 
to the flow associated with the "average supergranule" as seen by, e.g. |Duvall| 
and Hanasoge (2012). 
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Figure 1. Two-dimensional vector plots of the model supergranulation flow field. The top 
panel shows a horizontal slice through the horizontal flow field at the photospheric level 
(z = 0) and the bottom panel shows a vertical slice through the flow field at the center of the 
supcrgranule (y = 0). The full computational domain extends from -50 Mm to +50 Mm in the 
horizontal directions, where (x, y) = (0, 0) corresponds to the center of the supergranulation 
cell. Along the photospheric plane, the maximum horizontal and vertical velocities are 250 m 
s and 20 m s . 



From Equation ([I]), the radial and vertical mass-flux distributions, respec- 
tively, are 

P Vr = -^f, (4) 

We chose to model the v r component of flow as a sum of outflow [u ut] an< ^ 
inflow [v[ n ] Gaussian distributions, such that 

(z — zi) 2 (z — Z2) 2 
v r = %ut + v m = [«i CX P( £ja ) - Q 2 exp( )]F(r), (6) 

where a,\ — 250 m s _1 , a 2 — 6.27 m s^ 1 , Z\ — 0.2 Mm, z 2 = —15 Mm, and D\ = 
D 2 = 5 Mm. Using mass conservation, the coefficient a 2 is calculated from the 
known density distribution p and w ut with chosen coefficient at±. The potential 
function can then be calculated directly from 

Mr>*) = ~ f P(C)«r(C)dC, (7) 

J Zbot 
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In obtaining Equation |7|) , we set the integration constant equal to zero since we 
want the vertical velocity v z to be zero at the lower boundary (see Equation [8] 
below). 

With A$ fully prescribed, the vertical flow [v z ] is calculated from Equation 

As dAfh . , xr F(r) dFir), 

P^ = - + ^T=^W— + (8) 
r or r or 

Note that due to the large density gradient with solar depth, outflow velocities 
near the photosphere are much larger than return flow velocities at greater depth 
in order to conserve mass flux pv (Figure [T]) . 



3. Travel-time Measurements 



Starting from the vertical velocity taken at a height of 200 km, we used he- 



lioseismic holography ( Lindsey and Braun 2000 ) to measure center-quadrant 



local-control correlations. These correlations are analagous to center-quadrant 
time-distance correlations that have traditionally been used to measure time- 
distance travel times (e.g. Gizon and Birch 2005 Gizon, Birch, and Spruit 



2010). We used the quadrant geometry and ridge filters described by Braun and 



Birch (20081. In addition to the ridge filters, we also applied band-pass frequency 



filters to isolate particular ranges in frequency. The filters had central frequencies 
of 2.75, 3.00, 3.25, 3.50, 3.75, 4.00, 4.25, and 4.50 mHz and all had widths of 
0.25 mHz. We use Sr x (r) (§T y (r)) to denote x (y) direction center-quadrant 
travel-time shifts and (Sr ; = <5r out (r) — Sri n (r) to denote out minus in center- 
quadrant travel-time shifts, where 6T on t(v) (<5ri n (r)) represent travel-time shifts 
of the observed outgoing (incoming) waves. 

Figure [2] shows the power spectrum of the vertical velocity at a height of 
200 km above the photosphere in the simulation without any imposed flows. 
The resonant frequencies of the simulation are close to those from model S. As 
a result, we apply the ridge filters and holography Green's functions of |Braun] 
and Birch (20081 without modification. The travel-time measurements shown 
in Figures |3jand|4j were produced by subtracting the travel-time measurements 
from the simulation without flows (i.e., "noise subtraction"; Werne, Birch, and 



Julien 20041 



Figure |3| shows center-quadrant travel-time differences [5t x ] for a range of ra- 
dial orders and frequencies for the case where the synthetic data is produced from 
the supergranulation model. The flow produces travel-time shifts that generally 
decrease in amplitude with increasing radial order and increase in amplitude 
with increasing frequency (i.e., the travel-time shifts decrease with increasing 
horizontal phase speed). The spatial pattern is much the same in all cases and 
is what would be expected given the flow geometry in Figure [T] 

Figure [4] shows the <5r i travel-time differences resulting from the supergranu- 
lation model. The general trend of the travel-time shifts with radial order and fre- 
quency band is consistent with that in Figure |3j although the 6r i measurements 
feature a larger signal-to-noise ratio. 
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0.5 1 1.5 2 2.5 3 

k [rad/Mm] 

Figu re 2. Power spe ctrum for waves propagated through the CSM_A background model 
(Schu nker et aL]|2011| . The solid black lines denote the computed quiet Sun Model S mode 
frequencies and the dashed line denotes the horizontal phase speed corresponding to the sound 
speed at the bottom of the computational box. 



4. Flow Inversions 

The inverse problem is to estimate the flow field defined by the supergranule 
model given the travel-time measurements in Section [3j Here we show some 
example inversions in which we use the St x , Sr y , and St 01 maps to infer the 
subsurface flows v. 

We carry out the inversions of the travel-time maps with a noise level cor- 
responding to what we would expect for an average over 100 supergranules 
measured for 24 hours each. To generate the input travel-time maps we start 
from the noise-subtraction travel-time maps (Figures. [3] and [4} and add noise 
(computed from the simulations carried out in the flow-free reference model) 
with an amplitude reduced by 1/Vl00. 

We use kernels computed in the Born approximation using the method of 
Birch and Gizon (2007) with erratum |Birch, Gizon, and Burston (2011), and 



account for the holography Green's functions using the approach of |Birch et al. 
(2011). The vector- valued kernel functions K relate the flows in the interior to 
the travel-time maps, 

s n(r) = JJJ Ki(r' - r, z) ■ v(r', z) dr'dz . (9) 

In the above equation, we have separate maps and kernel functions for each 
combination of ridge filter and frequency range. 

Figure [5] compares the measured travel-time shifts with the travel-time shifts 
that we would expect from Equation ^ evaluated for the known flow in the 
simulation. At low frequencies there is general, although not perfect, agreement. 
At high frequencies, especially 4.5 mHz, there are very significant differences. 
This disagreement may be a result of the simplifications made in the calculation 
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Figure 3. Noise subtracted center-quadrant travel-time differences 5r x (r) calculated using 
ridge filters and 0.25 mHz wide band-pass frequency filters for the supergranulation simulation. 
The label at the bottom of each column denotes the radial order of the filter {e.g. n = is the 
/-mode, n = 1 is the pi-mode, etc.) and the factor by which the mode data has been scaled 
(e.g. x3 means the data has been multiplied by a factor of three). 



of the sensitivity kernels, for example the neglect of the sponge layers in the 
simulation. We discuss this issue in more detail in Sections <|5]and[6] 
In order to discretize Equation ([9]), we parameterized the flow [v] as 

N 

v(r,z) = j2Y, a ^(zy k - r (10) 

3=1 k 
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Figure 4. Noise subtracted center-annulus travel-time differences <5r j calculated using ridge 
filters and 0.25 mHz wide band-pass frequency filters for the supergranulation simulation. The 
label at the bottom of each column denotes the radial order of the filter (e.g. n = is the 
/-mode, n = 1 is the pi-mode, etc.) and the factor by which the mode data has been scaled 
(e.g. x3 means the data has been multiplied by a factor of three). The signal-to-noise ratio of 
the travel-time differences are generally greater than those of St x (Figure [3j. 



where the sum over j is over a total of N basis functions 4>j( z ) an d the sum 
over k is a two-dimensional inverse FFT on the grid that the travel-time maps 
are measured on (the simulation domain is periodic; as a result, the travel-time 
maps are periodic as well). Here we choose the basis functions <t>j{z) to be one 
when Zj < z < Zj+i and zero otherwise. The points Zj consist of a grid of 56 
points that are equally spaced in acoustic depth and cover the range from 0.5 
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Figure 6. Demonstration of different regularization applied to the v x inversion with corre- 
sponding decrease in resolved structure. The left panel shows the true v x at z = —0.8 Mm. 
From left to right, the remaining panels show the inversion results for v x for regularization 
parameters (A?,A|/fe^) = (1.0 X 10~ 8 , 5.66 X 10~ 5 ) Mm -1 (ms -1 ) -2 (second panel from 
left), (A 2 ,A|/h|) = (3.24 X 10~ 8 ,1.76 X 10~ 4 ) Mm" 1 (m s" 1 )" 2 (second panel from right), 
(A 2 , A 2 //i 2 ) = (1.76 X 10~ 7 , 9.31 Xl0" 4 ) Mm" 1 (m s" 1 )" 2 (right panel), where h x = 1/3 Mm 
is the grid spacing. The upper (lower) bound on the color scale in each panel is 150 (-50) m s — , 
with black (white) corresponding to positive (negative) values. 



Mm above the photosphere to 10 Mm below the photosphere. The grid spacing 
in fc-space is hk — 0.063 rad Mm -1 . 

We used the MCD approach (iJacobsen et al. 1999 1 to split the full inversion 



problem into small one-dimensional (z-only) inversion problems at each horizon- 
tal wavenumber. These problems we solved using RLS (e.g. Zhao, Kosovichev, 
and Duvall 2001 Kosovichev 1996 in the context of local helioseismology) with 



a regularization term given by the depth integral (at each horizontal wavevector 
fc) of the quantity v 2 x + v 2 + 10u|; the factor ten was chosen to reflect the 
preconception that vertical velocities are on average smaller than horizontal 
flows. Following Couvidat et oT](2005|, we used a fc-dependent regularization 
parameter with A 3 " 



\k\\ 2 , where (Ai,A2) are fixed parameters. As the 
regularization parameters are increased, the amplitude of the large-scale flow in 
the inversion result is reduced. We chose the parameters by inspection for each 
inversion based on a compromise between the desire to average out small-scale 
noise and the goal of retaining the large-scale flow features. Figure [6] shows a 
comparison of different regularization parameters applied to an inversion result 
for a portion of the v x velocity field. 

The inversion results presented herein reflect an estimate of the mean flow 
averaged over 100 supergranules for a 24-hour window. Application to real data 
necessitates consideration of the effective averaging window relative to the time 
scale over which the coherent structure evolves; a typical supergranule lifetime 



2008). 



is estimated to be a day or more (Worden and Simon 1976 Hirzberger et al. 



5. Results 

Figures [7] and [8] show inversion results using all available n = 0, n = 1, and n = 2 
measurements for (Af , A|//i|) = (3.24x 10" 8 , 1.69 x 10~ 4 ) Mm" 1 (ms^ 1 )- 2 com- 
pared to the true v x and v z velocity components in the supergranule simulation, 
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respectively. The structure of the i^-inverted fields in Figure [7] shares qualitative 
agreement with the true flow, although misses the signal at the top and bottom 
(below about z = —4 Mm) of the physical domain. The underestimation of 
flow magnitude in the lower layers is consistent with prior work (e.g. Zhao, 



Kosovichev, and Duvall 2001 ). The v z inversion in Figure [8] fails to qualitatively 
capture the true flow throughout the physical domain. As we will discuss later, 
this is due to cross-talk effects between components of the velocity. 
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Figure 7. Result of hclioscismic inversion for the x-component of the flow in the supergranulc 
simulation using all available n = 0, n = 1, and n = 2 measurements. The top-left plot 
compares the true (top panel) and inverted (bottom panel) v x (x,y) flow at z = —0.8 Mm. 
The bottom- left plot compares the true (left panel) and inverted (right panel) v x {x, z) flow 
at y = Mm. The plots on the right compare the centerline true (solid line) and inverted 
(dashed line with error bars) ^-velocity at z = —0.8 Mm (top panel), z = —2.0 Mm (middle 
panel), and z = —3.2 Mm (bottom panel). The horizontal axes are cropped relative to the full 
computational domain. The inversion produces velocity magnitudes that are comparable to 
the true signal, however misses the velocity signal at the top and lower portion of the domain. 



We tested the importance of including the n = 2 measurements by recom- 
puting the inversions shown in Figures [7] and [8] using only the n = and 
n = 1 measurements. There was no discernible visual difference in the inversion 
results for the horizontal and vertical velocity due to exclusion of the n = 2 
measurements, which indicates that the RLS technique is relatively insensitive 
to the addition of higher radial order measurements. 

We tested the effect of noise in the measurements by recalculating the flow 
inversions using only the n = and n = 1 measurements shown in Figures [3] 
and |4j which excludes the noisiest measurements (i.e. the whited-out boxes in 
the grid in Figures [3] and [4} that were included in the inversion calculations 
shown in Figures [7] & [5J There was no discernible visual difference in the in- 
version results for the horizontal and vertical velocity due to exclusion of the 
noisiest measurements, which supports the conclusion that the RLS technique 
is relatively insensitive to changes in the number and quality of measurements 
used in the inversion calculations. 
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Figure 8. Result of helioseismic inversion for the z-component of the flow in the supergranule 
simulation using all available n = 0, n = 1, and n = 2 measurements. The top left plot 
compares the true (top panel) and inverted (bottom panel) v z (x,y) flow at z = —0.8 Mm. 
The bottom left plot compares the true (left panel) and inverted (right panel) v z (x, z) flow 
at y = Mm. The plots on the right compare the centerline true (solid line) and inverted 
(dashed line with error bars) ti z -velocity at z = —0.8 Mm (top panel), z = —2.0 Mm (middle 
panel), and z = —3.2 Mm (bottom panel). The horizontal axes are cropped relative to the full 
computational domain. The inversion fails to qualitatively capture the true v z flow throughout 
the physical domain. As we will discuss later, this is due to cross-talk effects between signals. 



We calculated center-quadrant travel-time differences from Equation (J9j) using 
the kernels and the inferred flow from Figures [7[ and [8] and compared these to 
the measurements from the simulation. Figure p\ shows these comparisons and 
the calculated residual (difference between the two) for every combination of 
frequency filter, radial order, and measurement geometry. There is generally good 
correspondence between the perturbations in the n — and n = 1 measurements; 
however, some of the 5t x and 5r y measurements, as well as all of the n = 2 
measurements, contain significantly more noise. In some of the higher-frequency 
measurements, there is noticeable structure in the residual. This may be due 
to inaccuracies in the kernels (see Figure [5]) . We tested the sensitivity of the 
inversion to the high-frequency measurements by repeating the inversion using 
the same regularization parameters, but with the measurements at 4.25 mHz 
and 4.5 mHz removed. The main impact of removing these measurements was a 
reduction in the amplitude of the inferred vertical flow by about a factor of four. 

Our findings are generally consistent with those of Zhao et al. ( 2007 1 , who used 
the time-distance method and kernels computed from the ray approximation to 
perform inversions of simulated supergranule-scale convection. They were able to 
obtain reasonable inversion results only down to about 4 Mm below the surface, 
although they credit this to a limit in the largest annulus radius used. Here, 
the sensitivity functions for many of the travel-time measurements extend well 
below z = —4 Mm. This suggests that, with a sufficiently small noise level, 
it will be possible to detect flows below this depth (although with a spatial 
resolution that decreases with increasing depth) . The structure of the horizontal 
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Figure 10. Horizontally integrated 5t x kernels grouped by radial order and frequency, nor- 
malized by their respective rms noise level. Shown are n = kernels (solid line), n = 1 kernels 
(dashed lines), and n = 2 kernels (dot-dashed lines). The left- most panel shows a slice through 
the v x inversion result for reference. The figure demonstrates that the lobed nature of the v x 
inversion is due to the form of the kernels, which exhibit localized minima and maxima as a 
function of z. Furthermore, the sensitivity of the n = 2 kernels is generally less than the n = 
and n = 1 kernels, which means that use of higher radial-order measurements in the inversion 
calculations will have decreasing effectiveness on the result. 



velocity inversion results can be better explained by observing the form of the 
sensitivity kernels. Figure 10 contains a vertical slice through the v x inversion 
(left-most panel) and horizontally integrated St x kernels grouped by radial order 
and frequency, and normalized by the rms noise level. The n = kernels have a 
single band of sensitivity with a maximum located in between and -1 Mm. The 
n = 1 kernels have two regions of sensitivity which explains why the v x inversion 
has two local maxima in the vertical direction. The n = 2 kernels have several 
regions of sensitivity; however, inclusion of the n — 2 measurements did not 
significantly effect the performance of the inversions because the signal-to-noise 
ratio of the n = 2 measurements is about a factor of two lower than that of the 
n = and n = 1 measurements (Figure |4| . 

Figure [TT] shows the averaging kernel that relates the inferred v x to the true 
.T-componcnt of the velocity. The RLS inversion procedure is able to produce rea- 
sonably localized ^-component averaging kernels with minimal leakage outside 
the target region. 



Zhao et al. ( 2007 1 discuss the influence of cross talk on the inversions for the 



vertical velocity. Specifically, they describe a region just below the photosphere 
where the v z inversion is of the incorrect sign. We find a similar phenomenon 
in our results, illustrated by the band of negative vertical velocity near z = 
Mm (Figure |8]). Here we will use the notation v t j to denote the contribution of 
the j component of the (known) model velocity to the inversion result for the 
i component of the velocity; computing each of the Wy involves the horizontal 
convolution of the appropriate averaging kernel with a particular component 
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Figure 11. Averaging kernel for the v x inversion for z = —2.0 Mm. The left plot shows 
the two-dimensional slice through the kernel at z = —2.0 Mm and the right plot shows the 
horizontally-integrated vertical distribution of the kernel. 



of the known model velocity (see, e.g. Jackiewicz, Gizon, and Birch 2008 for 



examples). Figure 12 shows the contributions of each of the components of the 
model flow to the inferred vertical flow (i.e. v zx , v zy , v zz ). The contribution of 
the x- and y-components of the true flow to the inferred vertical flow are as 
large as the contribution from the true vertical flow. The sum of these velocities 
(vtot = v zx + v zy + v zz ) is qualitatively similar to the vertical flow inferred from 
the inversion. If the kernels were exactly correct and there was no noise then v to t 
would be identical to the inferred vertical flow. This computation highlights the 
importance of controlling the cross -talk in inversions for velocity (this has been 



achieved using OLA inversions by Jackiewicz, Gizon, and Birch 



2008 



Svanda 



et al. 2011). 



6. Discussion 



We computed RLS inversions from travel-time measurements of simulated data 
created by numerically propagating waves through a simple supergranule-like 
flow field. The benefit of this approach is the opportunity for comparison with a 
known flow field for direct evaluation of the performance of the measurement and 
inversion techniques. Furthermore, the use of simulated data allows for isolation 
of complicating effects. For instance, here travel-time differences are due solely 
to the supergranule flow; travel-time differences calculated from observational 
solar data contain the effects of many more physical variables. 

We found that the inversion of the v x - and Wj,-velocity shared general features 
of the true flow reasonably well down to a depth of about 4 Mm below the 
surface. The inversion of the vertical velocity performed poorly throughout the 
domain, particularly near the surface where the computations produced a result 
that was of incorrect sign. 

Underestimation of large-scale signal strength may be attributed to choice of 
regularization, which averages out noise but also causes some loss of flow ampli- 
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x[Mm] 

Figure 12. Convolutions between the ^-directional kernel and components of the model 
flow [v zx , squares; v zy , circles; v zz , triangles], the summation of the convolved components 
[vtot = v zx + v zy + v zz , dot-dashed line], the inversion result [v z inv, dashed line], and the 
model supergranulc flow [v z sim, solid line]. The v zx - and -u zy -components have a large effect 
on the calculation of the total signal [vtot] , which demonstrates that much of the error in the 
u z -inversion can be attributed to cross-talk effects. 

tude. We tested the inversions over a broad range of rcgularization parameters, 
and the results presented here reflect the best scenarios. 

We found the inversions, and hence the averaging kernels, to be relatively 
irresponsive to an increase in the number of input travel-time maps or to ex- 
clusion of noisy travel-time maps. This is a credit to the capability of the RLS 
technique in minimizing the effects of noise; however, the significant error in the 
inversion results suggests that the forward modeling is inadequate. 

The sensitivity kernels relate variations in travel-time shifts to variations in 
the solar model; capturing this sensitivity accurately is dependent on making 
appropriate assumptions regarding the physical state (in our case, the simulated 
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state) and the measurement procedure (e.g. Gizon and Birch 2002 Birch, Koso 



vichev, and Duvall 2004). As a test of the kernels, we convolved the resulting 
averaging kernels with the components of the true supergranule flow field, and 
compared these results with the respective inversion results. We found that the 
convolutions contained qualitatively the same structure (and hence comparable 
errors) as the inversion results. 

There is very little sensitivity to flows above the photosphere because only 
waves of frequency comparable to the acoustic cutoff frequency (about 5.5 mHz) 
sample this region. The calculation of the kernels uses a zero Lagrangian pres- 
sure perturbation upper boundary condition, an assumption which is likely not 
accurate for high-frequency waves. Furthermore, the simulation uses absorbing 
sponges at the top and bottom boundaries, which may affect the wave dynamics 
significantly differently from what is accounted for in the kernel calculations. 
Perhaps for these reasons, the inversions are unable to reproduce the flow struc- 
ture near the top of the computational domain. Other constraints on the forward 
modeling and approximations used in the design of the sensitivity kernels may 
need to be considered in order to improve the performance of flow inversions. 

In conclusion, our findings indicate that errors in the RLS inversions result 
from the combined effects of cross talk between signals, choice of regularization 
parameters, design of sensitivity kernels, and inconsistencies among boundary 
conditions, all of which are open areas of research requiring further study. We 
note that improvement in performance, particularly in regard to cross-talk ef- 
fects, has been seen in the optimally localized average (OLA) approach to reduce 



Jackiewicz, Gizon, and Birch 


2008 


Svanda et al. 


2011 


Jackiewicz et al. 2012 



Acknowledgements DED and ACB acknowledge support from NASA contracts NNH09CF68C| 
and NNH07CD25C. DCB acknowledges support from NASA contract NNH09CE41C. The 
data from the simulations used in this study is located on a publicly accessible website 
(http://www.cora.nwra.com/~dombroski/supergranulation/) to facilitate further investigation. 
The authors thank the referee for helpful suggestions. 



References 

Birch, A.C., Gizon, L.: 2007, Linear sensitivity of helioseismic travel times to 



local flows. Astronom. Nach. 328(3-4), 228-233. doi 10. 1002/asna. 200610724. 
http://dx.doi.org/10.1002/asna.200610724. 

Birch, A.C., Gizon, L., Burston, R.: 2011, Erratum: Linear sensitivity of 
helioseismic travel times to local flows. Astron. Nach. 332(6), 658-658. 
doi |10.1002/asna.201111557i http://dx.doi.org/10.1002/asna.201111557. 

Birch, A.C., Kosovichev, A.G., Duvall, T.L. Jr.: 2004, Sensitivity of Acous- 
tic Wave Travel Times to Sound-Speed Perturbations in the Solar Interior. 



Astrophys. J. 608, 580-600. doi: 10.1086/386361 



SOLA: manuscript_arxiv.tex; 30 November 2012; 2:21; p. 18 



Testing Helioseismic Holography Inversions 



19 



Birch, A.C., Parchevsky, K.V., Braun, D.C., Kosovichev, A.G.: 2011, "Hare 
and Hounds" Tests of Helioseismic Holography. Solar Phys. 272, 11-28. 
ADS |2011SoPh..272...11B [ doi |10.1007/sll207-011-9799TLl 



Braun, D.C., Birch, A.C.: 2008, Surface-Focused Seismic Holography of 
Sunspots: I . Observations. Solar Phy s. 251, 267-289. ADS |2008SoPh..251.| 
|.267B[ doi |10.1007/sll207-008-9152-5| 

Braun, D.C., Birch, A.C., Lindsey, C: 2004, Local Helioseismology of Near- 
Surface Flows. In: Dancsy, D. (cd.) SOHO 14 Helio- and Asteroseismology: 
Towards a Golden Future, ESA SP-559, ESA, Noordwijk 337. 

Braun, D.C., Birch, A.C., Benson, D., Stein, R.F., Nordhmd, A.: 2007, He- 
lioseismic Holography of Simulated Solar Convection and Prospects for the 
Detection of Small-Scale Subsurface Flows. Astrophys. J. 669, 1395-1405. 
doi |10.1086/521782| 

Christensen-Dalsgaard, J., Dappen, W., Ajukov, S.V., Anderson, E.R., Antia, 
H.M., Basu, S., Baturin, V.A., Berthomieu, G., Chaboyer, B., Chitre, S.M., 
Cox, A.N., Demarque, P., Donatowicz, J., Dziembowski, W.A., Gabriel, M., 
Gough, D.O., Guenther, D.B., Guzik, J.A., Harvey, J.W., Hill, F., Houdek, 
G., Iglesias, C.A., Kosovichev, A.G., Leibacher, J.W., Morel, P., Profhtt, 
C.R., Provost, J., Reiter, J., Rhodes, E.J. Jr., Rogers, F.J., Roxburgh, I.W., 
Thompson, M.J., Ulrich, R. K.: 1996, The Current State of Solar Modeling. 
Science 272, 1286-1292. doi |10.1126/science.272.5266.1286[ 

Couvidat, S., Gizon, L., Birch, A.C., Larsen, R.M., Kosovichev, A.G.: 2005, 
Time-Distance Helioseismology: Inversion of Noisy Correlated Data. Astro- 



phys. J. 158, 217-229. doi: 10.1086/430423 



Crouch, A.D., Charbonneau, P., Thibault, K.: 2007, Supergranulation as an 
Emergent Length Scale. Astrophys. J. 662, 715-729. doi: 10.1086/515564. 

De Rosa, M.L., Toomre, J.: 2004, Evolution of Solar Supergranulation. Astro- 
phys. J. 616, 1242-1260. doi: 10.1086/424920| 

Duvall, T.L. Jr.: 1980, The equatorial rotation rate of the supergranula- 
tion cells. So lar Phys. 66, 213-221. ADS |1980SoPh 66..213D1 doi |10.1007/| 
BF00150578. 



Duvall, T.L. Jr.: 1998, Recent Results and Theoretical Advances in Local Helio- 
seismology. In: S. Korzcnnik (cd.) Structure and Dynamics of the Interior of 
the Sun and Sun-like Stars, ESA SP-418, ESA Noordwijk 581. 

Duvall, T.L. Jr., Birch, A.C.: 2010, The Vertical Component of the Supergranular 
Motion. Astrophys. J. 725, L47-L51. doi |10.1088/2041-8205/725/l/L47| 

Duvall, T.L., Hanasoge, S.M.: 2012, Subsurface Supergranular Vertical Flows as 
Measured Using Large Distance Separations in Time-Distance Helioseismol- 
ogy. Solar Phys., 136. doi |10.1007/sll207-012-0010-0| 



SOLA: manuscript_arxiv.tex; 30 November 2012; 2:21; p. 19 



20 



D.E. Dombroski et al. 



Duvall, T.L. Jr., Jefferies, S.M., Harvey, J.W., Pomerantz, M.A.: 1993, T ime- 
distance helioseismology. Nature 362, 430-432. doi |10.1038/362430a0[ 

Gizon, L., Birch, A.C.: 2002, Time-Distance Helioseismology: The Forward 
Problem for Random Distributed Sources. Astrophys. J. 571, 966-986. 
doi |10.1086/340015| 

Gizon, L., Birch, A.C.: 2005, Local Helioseismology. Liv. Rev. in Sol. Phys. 2, 6. 
url: http://solarphysics.livingreviews.org/Articles/lrsp-2005-6/ (cited in October 
2011). 

Gizon, L., Birch, A.C., Spruit, H.C.: 2010, Local Helioseismology: Three Dimen- 
sional Imaging of the Solar Interior. Annu. Rev. Astron. Astrophys. 48, 289. 
doi |10.1146/annurev-astro-082708-101722| 

Hagenaar, H.J., Schrijver, C.J., Title, A.M.: 1997, The Distribution of Cell Sizes 
of the Solar Chromospheric Network. Astrophys. J. 481, 988. doi: 10.1086/ 
13040661 

Hanasoge, S.M., Duvall, T.L.: 2006, Solar acoustic simulator. In: Proceedings of 
SOHO 18/GONG 2006/HELAS I, Beyond the spherical Sun, ESA SP-624. 
ESA Noordwijk, 40. 

Hanasoge, S.M., Duvall, T.L. Jr., Couvidat, S.: 2007, Validation of Helioseismol- 
ogy through Forward Modeling: Realization Noise Subtraction and Kernels. 
Astrophys. J. 664, 1234-1243. doi |10.1086/519070| 

Hanasoge, S.M., Larsen, R.M., Duvall, T.L. Jr., De Rosa, M.L., Hurlburt, 
N.E., Schou, J., Roth, M., Christensen-Dalsgaard, J., Lele, S.K.: 2006, 
Computational Acoustics in Spherical Geometry: Steps toward Validating 
Helioseismology. Astrophys. J. 648, 1268-1275. doi |10. 1086/ 505927 

Hanasoge, S.M., Couvidat, S., Rajaguru, S.P., Birch, A.C.: 2008, Impact of 
locally suppressed wave sources on helioseismic traveltimes. Mon. Not. Roy. 
Astron. Soc. 391, 1931-1939. doi |10.1111/j.l365-2966.2008.14013.x| 

Hart, A.B.: 1954, Motions in the Sun at the photospheric level. IV. The equato- 
rial rotation and possible velocity fields in the photosphere. Mon. Not. Roy. 
Astron. Soc. 114, 17. 

Hart, A.B.: 1956, Motions in the Sun at the photospheric level. VI. Large-scale 
motions in the equatorial region. Mon. Not. Roy. Astron. Soc. 116, 38. 

Hathaway, D.H.: 1992, Spherical harmonic analysis of steady photospheric 
flows. II. So lar Phys. 137, 15-32. ADS |1992SoPh..l37...15H| doi |10.1007/| 
IBF001465731 

Hathaway, D.H., Beck, J.G., Han, S., Raymond, J.: 2002, Radial Flows in Super- 
granules. Solar Phys. 205, 25-38. ADS |2002SoPh 205...25H | doi |10.1023/Ai1 
101388121327S 



SOLA: manuscript_arxiv.tex; 30 November 2012; 2:21; p. 20 



Testing Helioseismic Holography Inversions 



21 



Hirzberger, J., Gizon, L., Solanki, S.K., Duvall, T.L.: 2008, Structure and 
Evolution of Supergranulation from Local Helioseismology. Solar Phys. 251, 
417-437. ADS:2008S oPh..251..417H] doi jl0.1007/sll207-008-920T8l 

Jackiewicz, J., Gizon, L., Birch, A.C.: 2008, High-Resolution Mapping of Flows 
in the Solar Interior: Fully Consistent OLA Inversion of Helioseismic Travel 



Times. Solar Phys. 251, 381-415. ADS:2008SoPh..251..381J doi: 10.1007/ 
|sll207-008-9158-zl 

Jackiewicz, J., Birch, A.C., Gizon, L., Hanasoge, S.M., Hohage, T., Rufho, J.B., 
Svanda, M.: 2012, Multichannel Three-Dimensional SOLA Inversion for Local 
Helioseismology. Solar P hys. 276, 19-33. ADS |2012SoPh..276...19J] doiJTU] 
|1007/sll207-011-9873-8| 

Jacobsen, B., Moller, I., Jensen, J., Efferso, F.: 1999, Multichannel deconvolu- 
tion, MCD, in geophysics and helioseismology. Phys. and Chem. of the Earth 
A 24, 215-220. doi |10.1016/S1464- 1895(99)0002T6| 

Kosovichev, A.G.: 1996, Tomographic Imaging of the Sun's Interior. Astrophys. 
J. 461, L55. doi: 10.1086/309989. 

Kosovichev, A.G., Duvall, T.L. Jr.: 1997, Acoustic tomography of solar con- 
vective flows and structures. In: Pijpers, F.J., Christensen-Dalsgaard, J., & 
Rosenthal, C.S. (ed.) SCORe'96 : Solar Convection and Oscillations and their 
Relationship, Astrophys. and Space Science Lib. 225, 241-260. 

Leighton, R.B., Noyes, R.W., Simon, G.W.: 1962, Velocity Fields in t he Solar 
Atmosphere. I. Preliminary Report. Astrophys. J. 135, 474. doi; 10.1086/ 
11472851 

Lele, S.K.: 1992, Compact Finite Difference Schemes with Spectral- like Resolu- 
tion. J. of Comp. Phys. 103, 16. doi |10.1016/002T"9 991(92)90324^R] 

Lindsey, C, Braun, D.C.: 1990, Helioseismic imaging of sunspots at their an- 



tipodes. Solar Phys. 126, 101-115. ADS:1990SoPh.. l26..101L[ doi{l0.1007/ 
IBF0015830T1 

Lindsey, C, Braun, D.C.: 2000, Basic Principles of Solar Acoustic Holography 
- (I nvited Review). Solar Ph ys. 192, 261-284. ADS:2000S oPh..l92" 261L] 
doi |10.1023/A:1005227200911| 

Ploner, S.R.O., Solanki, S.K., Gadun, A.S.: 2000, Is solar mesogranulation a 
surface phenomenon? Astron. Astrophys. 356, 1050-1054. 

Rast, M.P.: 2003, Supergranulation: new observation, possible explanation. In: 
Sawaya-Lacoste, H. (ed.) GONG+ 2002. Local and Global Helioseismology : 
the Present and Future, ESA SP-517, ESA Noordwijk, 163-172. 

Rieutord, M., Rincon, F.: 2010, The Sun's Supergranulation. Liv. Rev. So- 
lar Phys. 7, 2. http://solarphysics.livingreviews.org/Articles/lrsp-2010-2/ (cited 
October 2011). 



SOLA: manuscript_arxiv.tex; 30 November 2012; 2:21; p. 21 



22 



D.E. Dombroski et al. 



Rieutord, M., Roudier, T., Malherbe, J.M., Rincon, F.: 2000, On mesogran- 
ulation, network formation and supergranulation. Astron. Astrophys. 357, 
1063-1072. 

Rieutord, M., Roudier, T., Ludwig, H., Nordlund, A., Stein, R.: 2001, Are gran- 
ules good tracers of solar surface velocity fields? Astron. Astrophys. 377, 



L14-L17. doi 10.1051/0004-6361:20011160 



Rieutord, M., Roudier, T., Rincon, F., Malherbe, J., Meunier, N., Berger, T., 
Frank, Z.: 2010, On the power spectrum of solar surface flows. Astron. 
Astrophys. 512, A4. doi jlO. 1051/0004-6361/200913303| 

Schunker, H., Cameron, R.H., Gizon, L., Moradi, H.: 2011, Constructing and 
Characterising Solar Structure Models for Computational Helioseismology. 
Solar Phys., 124, 1-26. doi jl0.1007/sll207-011-9790^x 



Sekii, T., Kosovichev, A.G., Zhao, J., Tsuneta, S., Shibahashi, H., Berger, 
T.E., Ichimoto, K., Katsukawa, Y., Lites, B., Nagata, S., Shimizu, T., Shine, 
R.A., Suematsu, Y., Tarbell, T.D., Title, A.M.: 2007, Initial Helioseismic 
Observations by Hinode/SOT. Pub. Astron. Soc. Japan 59, 637-641. 

Simon, G.W., Leighton, R.B.: 1964, Velocity Fields in the Solar Atmosphere. 
III. Large-Scale Motions, the Chromospheric Network, and Magnetic Fields. 
Astrophys. J. 140, 1120-1149. doi |10.1086/148010| 

Svanda, M., Gizon, L., Hanasoge, S.M., Ustyugov, S.D.: 2011, Validated he- 
lioseismic inversions for 3D vector flows. Astron. Astrophys. 530, A148. 
doi jlO. 1051 /0004-6361 /201016426| 

Werne, J., Birch, A., Julien, K.: 2004, The Need for Control Experiments in Local 
Helioseismology. In: D. Danesy (ed.) SOHO 14 Helio- and Asteroseismology: 
Towards a Golden Future, ESA 5^-559, ESA Noordwijk, 172. 

Woodard, M.F.: 2007, Probing Supergranular Flow in the Solar Interior. 
Astrophys. J. 668, 1189-1195. doi:10.1086 /521391| 

Worden, S.P., Simon, G.W.: 1976, A study of supergranulation using a diode ar- 
ray magnetograph . Solar Phys. 46, 73-91. ADS |1976SoPh...46...73Wl doijlO] 



1007/BF00157555 



Zhao, J., Kosovichev, A.G.: 2003, On the inference of supergranular flows by 
time-distance helioseismology. In: H. Sawaya-Lacoste (ed.) GONG+ 2002. 
Local and Global Helioseismology: the Present and Future, ESA SP-517, ESA 
Noordwijk, 417-420. 

Zhao, J., Kosovichev, A.G., Duvall, T.L. Jr.: 2001, Investigation of Mass Flows 
beneath a Sunspot by Time-Distance Helioseismology. Astrophys. J. 557, 



384-388. doi: 10.1086/321491 



Zhao, J., Georgobiani, D., Kosovichev, A.G., Benson, D., Stein, R.F., Nordlund, 
A.: 2007, Validation of Time-Distance Helioseismology by Use of Realistic 
Simulations of Solar Convection. Astrophys. J. 659, 848-857. doi: 10.1086/ 
512009. 



SOLA: manuscript_arxiv.tex; 30 November 2012; 2:21; p. 22 



